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Abstract — In the literature, sealing laws for wireless mobile 
networks have been characterized under various models of 
node mobility and several assumptions on how communication 
occurs between nodes. To improve the realism in the analysis 
of scaling laws, we propose a new analytical framework. The 
framework is the first to consider a Levy flight mobility pattern, 
which is known to closely mimic human mobility patterns. 
Also, this is the first work that allows nodes to communicate 
while being mobile. Under this framework, delays (£>) to obtain 
various levels of per-node throughput (A) for Levy flight are 
suggested as D{\) — 0{^Ymn{n^+°'\, nP-)), where Levy flight 
is a random walk of a power-law flight distribution with an 
exponent a G (0, 2]. The same framework presents a new tighter 
tradeoff D{\) — 0(\/max(l, n\^)) for i.i.d. mobility, whose 
delays are lower than existing results for the same levels of per- 
node throughput. 

I. Introduction 

Since the work in [[T| that showed that mobility can be 
exploited to improve network throughput, there has been a 
plethora of work on this subject. A major effort in this 
direction has been in the design of delay tolerant networks 
(DTNs). However, this benefit in throughput comes at a sig- 
nificant delay cost. The amount of delays required to achieve a 
level of throughput for various mobility models such as i.i.d. 
mobility, random waypoint (RWP), random direction (RD), 
and Brownian motion (BM) have been extensively studied 
in ||2]-||6l. Specifically, the delay required for constant per- 
node throughput has been shown to grow as 8(n), which 
scales as fast as the network size n, for most mobility models 
including i.i.d. mobility, RWP, RD, and BM E), lE], Q, 
Despite significant advances in the development of delay- 
capacity scaling laws, there has been considerable skepticism 
regarding the applicability of the results to real mobile net- 
works because of various simplifying assumptions used in the 
analysis. 

In this paper, we address two issues towards making the 
delay-capacity tradeoff analysis more realistic: 1) contacts 
among nodes in the middle of their movements and 2) Levy 
mobility patterns of nodes in the network. In the literature, 
for mathematical simplicity, existing results have assumed that 

This work was supported in part by NSF grants CNS- 1065 136, CNS- 
1012700, CNS-0910868 and CNS-1016216, and ARO MURI Award 
W911NF-08-1-0238. 
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(a) Contacts while being mobile 



(b) A trajectory of Levy flight 



Fig. 1. Our considerations for improving the realism of delay-capacity 
tradeoff: (a) Conventionally, nodes A and B are assumed to have no contact 
opportunity, but we consider they can meet each other even during their 
movements, (b) We also consider that nodes follow Levy flight mobility. 



nodes show slotted movements, and they do not communicate 
with each other while being mobile. Thus, they make contacts 
with other nodes and transfer data only at the edge of time 
slots. In other words, as shown in Fig. [T] (a), the opportunity 
for meeting other nodes during mobility has been ignored, 
although such opportunities can substantially reduce packet 
delivery delays. Also, in this work we focus on the Levy 
flight model, which is widely accepted to closely mimic 
the actual movement of humans ||9l, ifTOl . The trajectory 
of this model is illustrated in Fig. [T| (b). To enhance the 
realism in the analysis of delay-capacity tradeoff, we develop 
a new analytical framework which takes both of these factors 
into account by developing a technique that characterizes the 
distribution of "first meeting time" among nodes conforming 
to Levy flight mobility in a two-dimensional space. It is 
important to note that the exact distribution of the first meeting 
time of Levy flight even in a one-dimensional space has 
been an open problem even though it has applicability in 
a diverse set of research problems (e.g., characterization of 
particle movements and animal movements) in physics and 
mathematics. It is also informative to note that the distribution 
of the first meeting time of BM, which can be considered as 
an extreme case of Levy flight, is also an open problem as 
noted in ifTH. ifTH. 

Levy flight, the mobility model we focus on in this paper, is 
a subset of Levy mobility in which a node moves from position 
to position in a constant time. Another special case. Levy 
walk, in which a node moves from one position to another 
in time proportional to the distance between the positions. 
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Except for the notion of the time required for each movement. 
Levy flight and Levy walk are fundamentafly the same random 
walk whose flight length distribution asymptotically follows 
a power-law fa{z) w where z and a denote the 

flight length (i.e., moving distance of each slotted movement) 
and the power-law slope ranging < a < 2, respectively. 
The heavy-tailed movements of Levy mobility render the 
delay characterization extremely challenging. Our framework 
addresses these challenges using theories from stochastic ge- 
ometry and probability, and provides a delay-capacity tradeoff 
for Levy flight. Also, for a simpler i.i.d. mobility model, we 
provide a tighter delay-capacity tradeoff compared to existing 
studies using the same framework. 

Fig.|2]and Table U summarize the new tradeoffs identified us- 
ing our analytical framework. The results show that the trade- 
off for Levy flight follows D{X) = 0(^min{ni+"A, n^}) to 
obtain a per-node throughput of A = 8(n~'') {0 < rj < 1/2) 
as shown in Figs. |3](a) and (b). These results are well aligned 
with the critical delay (i.e., minimum delay required to achieve 
A = uj{l/y/n)) suggested in ifTSl . Our tradeoffs show an 
important finding that the delay required to obtain constant 
per-node throughput (i.e., A = 8(1)) can be smaller than 
Q{n) in mobile networks with mobility models such as Levy 
flight with a < 1 and i.i.d. mobility. This is an important 
observation given that most of the existing studies present 
the delay required to obtain constant per-node throughput to 
be 8(?i) for almost all mobility models including the i.i.d. 
mobility. 

Our tradeoff for Levy flight becomes especially more inter- 
esting when we input a values from measurements presented 
in Table HI] into the tradeoff. This gives us a hint on how 
the performance of the network will scale in reality when 
the network consists of devices mainly carried (or driven) by 
humans. For a values between 0.53 and 1.81, the delays to 
obtain A = 8(1) are expected to lie between ©(n^ *^^) and 
0{n). This implies that in reality, a DTN mainly operated by 
human mobility may indeed experience less than 8(?i) delay 
in some areas. This observation of smaller delay suggests that 
mobile networks relying on opportunistic transmissions may 
have higher practical values in reality given that the delays 
have been overestimated by mobility and contact models with 
less realism. 

TABLE I 

The new delay-capacity tradeoffs identified using our 
analytical framework for levy flight and u.d. mobility 



D = 0(v/miii{rii+"A,n2}) D 
Q = 2 



Mobility 


Tradeoff D{X) 


_D(A) for A = 0(1) 


Levy flight 


0(Vmin(ni+"A,n2)) 


0(Vmin(ni+°=,n2)) 


i.i.d. 


0(Vmax(l,nA3)) 


0(nV2) 



The rest of the paper is organized as follows. We overview 
a list of related work in Section and introduce our system 
models and definitions of performance metrics in Section [III] 




D = 0(\/nifix{l,«A-'}) 



(a) Levy flight 



(b) i.i.d. mobility 



Fig. 2. (a) The delay-capacity tradeoffs for Levy flight with various a 
obtained from our analytical framework and (b) the new tradeoff for i.i.d. 
mobility (solid line) shown together with the previously known tradeoff 
(dotted line). 
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(a) A = e(l/v^) 



1 2 

(b) A = 0(1) 



Fig. 3. The upper bounds on D{X) for Levy flight with a S (0, 2] to obtain 
(a) A = e(l/Vn) and (b) A = 0(1). The critical delay for Levy flight 
derived in (13| is depicted in (a). 

TABLE II 

EXPERIMENTAL a VALUES FROM FIVE SITES PRESENTED IN [9). THESE 
VALUES ARE FROM EXPERIMENTAL GPS TRACES FROM FIVE SITES WITH 
MORE THAN 1 00 PARTICIPANTS IN TOTAL. 



Site 


a 


Site 


a 


Site 


a 


KAIST 
NCSU 


0.53 
1.27 


New York 
Disney World 


1.62 
1.20 


State fair 


1.81 



We then provide the intuition on how our analytical framework 
evaluates delay-capacity tradeoff using the properties of first 
meeting time of random walks in Section |IVl Based on the 
understanding in Section |IV] we analyze the tradeoffs of 
Levy flight and i.i.d. mobility models in Sections fV\ and [Vll 
respectively. After briefly concluding our work in Section [Vlll 
we provide full details of proofs used in Sections |V] and [Vll 
through Appendices. 

n. RELATED Work 

In lfT4l . it was shown that the per-node throughput of 
random wireless networks with n static nodes scales as 
(9(1/ v^). The result was later enhanced to Q{l/y/n) by 
exercising individual power control ifTSl . Grossglauser and 
Tse im proved that constant per-node throughput is achievable 
by using mobility when nodes follow ergodic and stationary 
mobility models. This contradicted the conventional belief that 
node mobility negatively impacts network capacity due to 
interruptions in connectivity. 

Many foflow-up studies H, [U, E], Q, d-Qll have 
been devoted to characterize and exploit the delay-capacity 
tradeoffs. In particular, the delay required to obtain constant 
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per-node throughput has been studied under various mobility 
models ||2l-||6l. The key message is that the delay of 2-hop 
relaying proposed in U] is Q{n) for most mobility models 
such as i.i.d. mobility, RD, RWP, and BM. 

The delay-capacity tradeoff for per-node throughput A = 
e(n-'') (0 < 7? < 1/2) is first presented in Eg as D{X) = 
0{nX) for RWP model. In @, the authors identified that 
D{\) ^ 6(max(l, nA^)) holds for i.i.d. mobility. Later, g] 
showed that D{X) = Q{n) holds for BM irrespective of A. 

More realistic mobility models. Levy mobility models 
known to closely capture human movement patterns, were 
first analyzed in ||T3| for a special case of the tradeoff. Using 
spatio-temporal probability density functions, the critical delay 
defined by the minimum delay required to achieve larger 
throughput than Q{l/y/n) is identified for Levy flight as well 
as Levy walk. 

Existing results on delay-capacity tradeoffs for mobile 
networks have been built under the assumption that nodes 
are able to communicate with each other only at the edge 
of time slots for slotted movements. Also, there has been 
no framework which can fully understand the delay-capacity 
tradeoff for Levy mobility. In this paper, we develop an 
analytical framework which handles both of these issues and 
are able to use this framework to characterize the delay- 
capacity tradeoff. 

III. System Model 

A. Network Model 

We consider a wireless mobile network indexed by n S 
N = {1, 2, . . .}, where, in the nth network, n nodes move on 
a completely wrapped-around disc X> (c K.^) whose radius 
scales as Without loss of generality, we set the 

radius and the center of the disc X> as y/n and = (0, 0), 
respectively, i.e., T> ^ {x E M."^ \ \x\ < \/n}. We assume that 
the density of the network is fixed to 1 as n increases 1 We 
also assume that all nodes are homogeneous in that each node 
generates data with the same intensity to its own destination. 
The packet generation process at each node is independent of 
node mobility. The generated packets are assumed to have no 
expiration until their delivery and the size of each node's buffer 
is assumed to be unlimited. Each packet can be delivered by 
either direct one-hop transmission or over multiple hops using 
relay nodes. 

To model interference in wireless networks, we adopt the 
following protocol model as in JS), lfT4l . Let Xi{t) denote the 
location of node i (= 1, . . . , n) at time t (> 0). Let i(i j) (t) = 
\Xi{t)—Xj{t)\ denote the Euclidean distance between nodes i 
and j at time t. Under the protocol model, nodes transmit 

'in all notations, a bold font symbol is used to denote a value or a set of 
values in M?. 

-This model is often referred to as an extended network model. In another 
model, called a unit network model, the network area is fixed to 1 and the 
density increases as n while the spacing and velocity of nodes scale as 1/ 



packets successfully at a constant rate W bits/sec, if and only 
if the following is satisfied; for a transmitter i, a receiver j 
and every other node u i,j transmitting simultaneously, 

hu,j)it) > (l + A)L(,,,)(t), 
where A is some positive constant. In addition, the distance 
between nodes i and j at time t should satisfy i(i j) (t) < r, 
where r (> 0) denotes the maximum communication range. 
We assume the fluid packet model 16], which allows concurrent 
transmissions of node pairs (with the rate divided by the 
number of pairs) interfering each other We denote by 11 
the class of all feasible scheduling schemes conforming our 
descriptions. 

B. Mobility Model 

In this subsection, we mathematically describe the Levy 
flight model and the i.i.d. mobility model. At time t = 0, 
node i chooses its location uniformly on the disc 7? (i.e., 
Xi{0) ^ Uniform(X')), which is independent of the others 
Xj{0) for j 7^ i. We assume that time is divided into slots of 
unit length and is indexed by fc G N. At the beginning of the 
fcth slot (i.e., at time t ^ k — 1), node i chooses its next slotted 
location Xi{k) according to the associated mobility model. 
During the kt\\ slot (i.e., during time t £ {k — l,fc]), node i 
moves from Xi{k—1) to Xi{k) with a constant velocity. Thus, 
Xi{k~l + 6) for 6 e (0, 1) is determined by Xi{k - 1) and 
Xi{k) as follows: 

X,{k -l + 5) = {l- 5)X,{k - 1) + 5X,{k). (1) 

Levy Flight Model. At the beginning of the fcth slot (i.e., at 
time t = k—1), node i chooses flight angle and flight length, 
denoted by 9i{k){e (0,27r]) and Zi{k) {> 0), respectively. 
During the fcth slot (i.e., during time t £ {k — l,fc]), node i 
moves from Xi{k — 1) to the selected direction 9i{k) of the 
distance Zi{k). Thus, the location Xi{k) is determined as 

X,{k)=X,{k-l) + V,{k), (2) 

where 

V,{k) ^ (Z,(fc)cos0,(fc),Z,(fc)sin0,(fc)). (3) 
The flight angle 9i{k) and the flight length Zi{k) are inde- 
pendent of each other and also independent of the previous 
locations Xi{t) for the times t G [0,fc — 1] before they 
are generated. Hence, Vi{k) is independent of Xi{t) for all 

t e [0, fc- 1]. 

Each flight angle di{k) and flight length Zi{k) are inde- 
pendent and identically distributed across node index i and 
slot index fc. Let 9 and Z he a generic random variable 
for 9i{k) and Zi{k), respectively. Then, the flight angle 9 is 
uniformly distributed over (0, 27r], and the flight length Z is 
generated from a random variable Z* having the Levy a-stable 
distribution 1201 by the relation Z = \Z*\. The probability 
density function of Z* is given by 
1 



2tt 



(4) 



4 



where ipz*{t) — E[e'*'^ ] is the characteristic function of Z* 
and is given by (pz*{t) = e^''''*!". Here, \s\ > is a 
scale factor determining the width of the distribution, and 
a e (0, 2] is a distribution parameter that specifies the shape 
(i.e., heavytail-ness) of the distribution. The flight length Z 
for a S (0, 1) has infinite mean and variance, while Z for 
a e [1,2) has finite mean but infinite variance. For a = 2, the 
Levy a-stable distribution reduces to a Gaussian distribution 
with a mean of zero and variance of 2s^, for which the flight 
length Z has finite mean and variance. 

Due to the complex form of the distribution, the Levy a- 
stable distribution for a E (0, 2) is often treated as a power-law 
type of asymptotic form: 



\l+a ' 



(5) 



where we use the notation a{z) ^ b{z) for any two real func- 
tions a{z) and b{z) to denote Mmz^aoiaiz) /b{z)] = 1 1201 . 
The form (|5]l is known to closely approximate the tail part of 
the distribution in and a number of papers in mathematics 
and physics, e.g., ET\ . Il22l . analyze Levy mobility using the 
form (|5]l. For mathematical tractability, in our analysis we will 
also use the asymptotic form (|5]l. Specifically, we assume that 
there exist constants c (> 0) and Zth (> 0) such that 

P{Z > z} = — , for all z > Zth. (6) 

i.i.d. Mobility Model. At the beginning of the fcth slot (i.e., 
at time t ~ k — 1), node i chooses Xi{k) uniformly on the 
disc I?, which is independent of its previous locations Xi{t) 
for the times < e [0, fc — 1] as well as the others Xj{k) for 
j 7^ i. Thus, Xi{k) is independent and identically distributed 
across node index i and slot index fc. 



C. Contact Model 

In our contact model, nodes are allowed to meet while being 
mobile. Hence, for a time t* in a domain {t\t > 0}, we say 
that nodes i and j meet at time t* (or are in contact at time t*) 
if they satisfy 

In the widely adopted contact model where nodes are 
allowed to meet only at the end of their movements (i.e., at 
slot boundaries), a meeting event can occur for a time fc* in a 
domain {fc | fc G {0} U N} satisfying X(ij)(fc*) < r. We call 
this class of contact model slotted contact model throughout 
the paper 

Mobile nodes are exposed to more contact opportunities in 
our contact model compared to the slotted contact model. 

D. Performance Metrics 

The key performance metrics of our interest are per-node 
throughput and average delay as defined next: 



Definition 1 (Per-node throughput). Let A.^-i{t) be the total 
number of bits received at the destination node i up to time t 
under a scheduling scheme tt G 11. Let Ajr be the per-node 
throughput under tt. Then, 



liminfiy^. 

t->oo n ^ — ' t 
1=1 



Definition 2 (Average delay). Let D-^ a^v be the time taken for 
the vth packet generated from the source node i to arrive at 
its destination node under a scheduling scheme tt G 11. Let 
Dt, be the average delay under tt. Then, 

1 " 1 



lim 



i—l v—1 

In this paper, we focus on analyzing the scaling property of 
the smallest average delay achieving per-node throughput A. 
We call this minimum average delay optimal delay throughout 
this paper. We focus on the throughput only in the range 
from 8(l/^/n) to 6(1), since this range corresponds to the 
case where mobility can be used to improve the per-node 
throughput. 

Definition 3 (Optimal delay). Let D{X) be the optimal delay 
to achieve per-node throughput A. It is then given by 

D(X) = inf D^. 

{Treni A„=A} 

E. Throughput Achieving Scheme 

We now consider a scheme tt that can realize per-node 
throughput Afr scaling from Q{l/^/n) to 0(1). The scheme tt 
operates as follows: 

• When a packet is generated from a source node and the 
destination of the packet is within the communication 
range of the source node, the packet is transmitted to the 
destination node immediately. 

• Otherwise, the source node broadcasts the packet to all 
neighboring nodes within its communication range. Note 
that this broadcast is only performed by the source node 
when the packet is generated. 

• Any nodes carrying the packet can deliver the packet to 
the destination node when they are within the communi- 
cation range of the destination node. 

• When one of the packets (including the original packet 
and the duplicated ones) reaches the destination node, all 
others are not considered for delivery. 

By appropriately scaling r as a function of n, the scheme tt 
can achieve per-node throughput A^r ranging from 6(l/\/n) 
to 8(1), as shown in the following lemma. 

Lemma 1. Let the communication range r scale as Q{n^) 
(0 < /? < 1/4). Then, the per-node throughput A# under the 
scheme tt scales as Q{n^^^). 

Proof: If the network has been running for a long enough time, 
all nodes become to work as relay nodes and begin to have 
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packets for all other nodes. Therefore, for a network with n/a„ 
disjoint area where a„ = 6(r^), all areas with more than two 
nodes can always be activated. Let denote the probability 
of having more than two nodes in an area. We then have 

— -n—[l 

n / n \ n 



1-1 



In addition, the total network throughput becomes 6„n/a„ and 
accordingly the per-node throughput is A^^ = fe„/a„. Without 
loss of generality, we assume r = (0 < /3 < 1/4). Then, 
the per-node throughput is given by 

A- = „-2^ - (1 - 7i2/3-i)" _ (1 _ „2^-i)"-i . 

In the following, we will show that 
A* 



lim 



-2P 



= lim {l-(l-7i2^i)"-n2^ (1-^2^1)""'} (7a) 



(7b) 



_ ( 1 - 2cxp(-l), if^ = 0, 

'jl, if /3e (0,1/4]. 

Hence, the per-node throughput A# under the scheme tt scales 
as 6(n~2/5). Note that 

lim log(l - n^'^-i)" = lim n^-^ log(l - n^'^-i)"'"'^ 

n— )-oo n— >-oo 

= lim n^^log(exp(-l)) 

n— >-oo 

^f-i, if/3 = 0, 
1-00, if /?e (0,1/4], 



which gives 
lim (1 — n 



cxp(-l), if/3 = 0, 



[0, if /3e (0,1/4] 

Similarly, 

lim log(n2'3(l-n2^-i)"-i) 

n— fcxD 

= lim {2/31ogn+(n-l)n2/5-ilog(l-n2^-i)"'""'} 

n— >-oo 

= lim {2/31ogn + n^^log(cxp(-l))} 

n— >^cxD 

^f-l, if/3 = 0, 
\-c», if /3g (0,1/4], 
which gives 



(8) 



.2/3n_^2^)-i^n-i ^ |exp(-l), if/3 = 0, 

lo, if /3e (0,1/4]. 



Iimn2^(l-n^^-i)^ 



(9) 



By applying ([8]l and (|9]) to dTaj i, we have ( |7b] i. This complete 
the proof. ■ 

Let D-nifi) be the average delay under tt when r ~ O(n^). 
Lemma [T| implies that by setting /? = — log„(A/A), the 
scheme tt achieves the per-node throughput A# = A. Since 
the scheme tt is of the class 11, the order of D{X) can be 
obtained from Djr{(3) with the use of /3 = — log„(V^), i.e., 

D{X)^ inf <^*(-log„(VA)). (10) 



IV. Preliminaries 

In this section, we provide the key intuition to understand 
how our analytical framework utilizes the properties of first 
meeting time in the derivation of delay-capacity tradeoffs 
under the Levy flight and the i.i.d. mobility models. We then 
sketch the challenges residing in our framework and briefly 
describe our approach to address these challenges. 

A. Delay Analysis with First Meeting Time 

The first meeting time of two nodes moving in a two- 
dimensional space, which is directly connected to Dj^, is 
defined below: 

Definition 4 (First meeting time). For i ^ j, the first meeting 
time of nodes i and j, denoted by T(^i.j), is defined as 
^inf{t >0|L(,,,)(t) <r}. 

Since T(^i.j} is independent and identically distributed across 
pair index (i, j), we use T to denote a generic random variable 
for 

Let £'(s,d) be a random variable representing the time taken 
by a packet generated from a source node s to arrive at a 
destination node d. Since the packet generation process is 
independent of node mobility, we consider that each packet 
is generated at time t ~ Q without loss of generality. Then, 
the packet delay under the scheme tt, denoted by can 
be expressed in terms of the first meeting time as 

fo, ifrfel(s), 

[min(T(i^d);i eT(s)), ifd^I(.s), 

where I{s) = {i \ i)(0) < r} denotes a set of node indices 
that are within the communication range of the node s at time 
t = 0. Note that s G X(s) by definition. Hence, the following 
equation represents the scheme tt described in Section |llll 

^7r:(s,d) = min (Tf^r.d) ; « e 2:(.s) \{d}). 
From Definition 2, the average delay Dji can be obtained by 

= E[min(%d);7eI(.s)\{4)]. (11) 

B. Distribution of the First Meeting Time 

In order to evaluate (fTTl i. the distribution of the first meeting 
time P{T(i < r} is essential. To obtain the distribution, we 
start from defining the following: let I(i.j) (k) {i ^ j,k £ N) 
be a random variable indicating the occurrence of a meeting 
event between nodes i and j during the fcth slot (i.e., time 
t e {k- l,fc]), i.e., 

'o, if L(i_j-){t) > r for all t e {k - l,k], 
1, if i(ij)(t) < r for some t <E (k — 1, k]. 
For notational simplicity, throughout this paper, we omit 
in and j)(-), unless there is confusion. We then 

define a function H{k, Iq) for fc G N and Iq G (r, 2y/n\, which 



D 
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denotes the probability that nodes i and j are not in contact 
during the fcth slot, conditioned on the fact that the initial 
distance between the nodes was Iq and after that the nodes 
have not been in contact by time t = k — 1, i.e., 

fp{/(i) = o|L(o) = ;o}, iffc==i, 

H{k,lo)^ <P{Iik) = 0\I{k-l) = ... 

[ =/(l) = 0,L(0) = /o}, iffc = 2,3,---. 

(12) 

Note that Iq is upper bounded by 2y/n since the radius of the 
disc X> is set to y/n. 

We find that the distribution of the first meeting time T 
can be obtained from the function H{k, Iq) as shown in the 
following lemma. 

Lemma 2. For r = 0, 1, 2, . . ., the complementary cumulative 
distribution function ( CCDF) of the first meeting time T can 
be obtained by 

P{T>t} = J ^ [l[H{k,lo))dF^o){lo): (13) 

where ^l(o)(') denotes the cumulative distribution function 
(CDF) of L{0), and we use the convention ^ ^■ 

Proof: For r = 0, the event {T > 0} implies the event 
{i(0) > r}, and vice versa. Hence, we have 

P{T > 0} = P{L(0) >r}= 1 dFLio) (lo)- 

Jr+ 

For r = 1, 2, . . ., the CCDF P{r > r} can be obtained by 

p{r > t} 

= P{L(0)>r,/(l) = ... = /(r)=0} 

= / P{/(l) = ... = /(r) = 0\LiO)^lo}dFno)ilo) 

Jr+ 

\\H{kM)) dFLio}{lo), 
fe=i 

which completes the proof. ■ 

The identity P{T > 0} = P{£(0) > r} shown in 
Lemma 12] has the following imphcations for P{T > 0}: (i) It 
is determined by the spatial distribution of nodes at time t = 
(which is assumed to be uniform on the disc T>). Hence, it is 
invariant for both the Levy flight and the i.i.d. mobility models, 
(ii) It represents the probability that two arbitrary nodes are out 
of the communication range at time t = 0. Since P{T > 0} is 
frequently used throughput this paper, we define it as Pq and 
summarize its implications using the following lemma. 

Lemma 3. Suppose that at time i = 0, nodes are distributed 
uniformly on a disc of radius ^/n. Define 

Po^P{r>0}(=P{L(0) >r}). (14) 

Then, Pq is bounded by 

1 - - < < 1 - 
11 6n 

Proof: See Appendix A. ■ 



C. Technical Challenge and Approach 

In our framework, characterizing the function H{k, Iq) 
in ([T2I1 which appears in the expression for P^T" > t} 
in ( fT3] l is the key to analyze the optimal delay. The major 
technical challenge arises from tracking meeting events in the 
middle of a time slot. The meetings over time are heavily 
correlated irrespective of mobility models. The correlation can 
be understood as follows: let us consider two consecutive slots, 
say the fcth and the (fc + l)th slots, for ease of explanation. 
The occurrence of a meeting event during the fcth slot (resp. 
the (fc + l)th slot) is determined by the locations of nodes 
i and j at the slot boundaries, i.e., at times t = fc — l,fc 
(resp. at times t = fc, fc + 1). Hence, both /(fc) and /(fc + 1) 
depend on the values of Xi{k) and Xj{k), and accordingly 
the sequence {/(fc)}^^^ is correlated in our contact model. 
Due to the complexity involved in the correlation, deriving 
the exact form of //(fc, /q) appears to be mathematically 
intractable. To address this challenge, we take a detour to 
derive a bound on //(fc, Iq) using theories from stochastic 
geometry and probability. The detailed analysis of //(fc, Iq) for 
the Levy flight model and the i.i.d. mobility model is presented 
in Lemmas [5] and [TOl respectively, which allow us to reach the 
conclusions of this paper. 

V. Delay Analysis for the Levy Flight Model 

In this section, we analyze the optimal delay under the Levy 
flight model. We use the following four steps in our analysis: 

• In Step 1, the average delay under our Levy flight model 
is formulated explicitly using the distribution of the first 
meeting time T. 

• In Step 2, we derive a bound on the distribution of T 
by characterizing the function //(fc, Iq) under the Levy 
flight model. The difficulty of handling contacts while 
being mobile is addressed in this step. 

• In Step 3, we connect the result of Step 2 to the delay 
scaling under the Levy flight model. 

• In Step 4, we study the delay-capacity tradeoff by com- 
bining the capacity scaling in Lemma [T| and the delay 
scaling obtained in Step 3. 

Step 1 (Formulation of the average delay using the 
first meeting time distribution): As shown in ( fTTl i. = 
E[min(T(j i G I(.s) \ {d})]. Under the Levy flight model, 
for i € I(s) are heavily correlated since the next slotted 
location Xi{k + 1) depends on the current location Xi{k). 
Note that all the nodes i <E /(s) are in proximity of the 
node s, and thus min(T(j ,i);i £ I{s) \ {d}) is not easily 
tractable. Therefore, we use the following bound to describe 
/)# using T. 

^# <E[r(3,d)] =E[r]. (15) 

For the simpler i.i.d. mobility model, we are able to derive a 
tighter bound on min(r(,;.rf); i G I(s) \ {d}). We present the 
result in Step 1 of Section IVII 
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Fig. 4. The illustration of S{lo). 

Let [T] denote the smallest integer greater than or equal 
to T. Then, since E[r] < E[[r]] < E[r] + 1, the order of 
E[[T]] is the same as that of E[r], and E[[r]] is an upper 
bound on as shown in the following lemma. 

Lemma 4. The average delay Dj^ of the scheme tt under the 
Levy flight model is bounded by 

Di,<E[\TW, (16) 

where T is the generic random variable for the first meeting 
time T^i,j) defined in Definition^ The expectation ^[[r]] can 
be obtained from the distribution of T by 

oo 

E[m] = Y^p{T>T}. 

T = 

Proof: From we have 5^ < E[T]. Since T < [T], 

we have E[r] < E[[r]], which gives ( fT6] l. Since the random 
variable \T~\ takes on only nonnegative integer values, the 
expectation E[[r]] can be obtained by 



E[m] 



E 

f=i 



P{rri >r} = ^p{r>f-i}, 



where the second equality comes from the property that 
P{ [T] > f} = P{T > f - 1} for all f 1, 2, . . .. Replacing 
f — 1 with T gives the lemma. ■ 

Step 2 (Characterization of the first meeting time distribu- 
tion): In this step, we first analyze the characteristics of the 
function H{k, Iq) under the Levy flight model (See Lemma|5]). 
By exploiting the characteristics, we then derive a bound on 
the first meeting time distribution (See Lemma|6l). This bound 
enables us to derive a formula for the expectation E[[T]] used 
in Lemma |4] (See Lemma |7]). 

Lemma 5. Under the Levy flight model, the function H{k, Iq) 
in ( 1721) has the following characteristics: 
(i) Let AV be a generic random variable for Vi{k) — 
representing a flight differential between nodes i and jj^ Then, 
geometrically the function H{1,Iq) can be viewed as the prob- 
ability of the flight differential falling into a set S(Iq) (c M?) 
defined as follows. Let'Dri'iJ') (c M^) denote a disc of radius r 
centered at u e M^, i.e., T>r{u) = {a; G | |a; - it| < r}. 

'The existence of the generic random variable for Vi{k) — Vj{k) is proven 
in Lemma [T4l in Appendix B. 



Let (v,w) denote a line connecting two points v,w E R^. 
For a fixed Iq £ {r, 2y/n\, define a set S{lo) as 

S{lo) ^{xe M^l (O,a;)nX»,.((0,-?o)) = 0}. (17) 
An example of Silo) is shown in Fig. |4] The set S{Iq) has a 
connection with the function H{1,Iq) as follows: 

H{l,lo)=P{AV eS{lo)}. 

(ii) The function H{1,-) is nondecreasing. 

(Hi) From (ii), we have H(l,lo) < P for all Iq S (r, 2-y/n], 

where 

P = H{l,2y^). (18) 

(iv) For k ~ 2,3, . . ., each function H{k, Iq) is also bounded 
above by P. Thus, for all k N and Iq e (r, 2^/n\, we have 

H(k,lo) < P. 

(v) There exist constants ci,Cu (> 0), and n,i, G N such that 
for all n > n,;,, P is bounded above and below by 



2ci 

P <1 

TT 



1 



P > 1 



2^. 

2a/2+2 



sm 



-1 / 



1 



■ -if r 



sm 



(19a) 
(19b) 



TT \2^/n — r ^ 
Proof: See Appendix B. ■ 
Based on the formula for P{T > t} in Lemma |2] and (iv) 
in Lemma [5] we derive a bound on PjT > t} in terms of P 
as follows: for t = 1, 2, . . ., we have from (iv) in Lemma |5] 
that 



0<l[Hik,lQ) < {P)\ 



(20) 



k=l 



Since (|20] | holds for all Iq e (r, 2-/ri\, by integrating ( l20b over 
Iq e (r, 2^/n], we have 

'XlH{k,lQ)) di^i(o)(/o) 

k=l 

< / (P)^dFi(o)(Zo) = (P)"Po. (21) 

Jr+ 

By combining (l2Tl i and Lemma |2] we have 

P{T>t} <{P)^Po, forr = 1,2,.... (22) 

Since P{T > r} = for r = 0, the bound in ^ also holds 
for T = 0. The above result is summarized in Lemma |6] 

Lemma 6. Under the Levy fiight model, the CCDF of the first 
meeting time T is bounded by 

P{T>t}<{P)-Po, forr = 0,1,..., 

where P and Pq are defined in and ( 1741 ), respectively. 

Proof: Combining Lemma |2] and (iv) in Lemma |5] gives 
Lemma |6] The detailed derivation was described earlier 
in (EOl-dlali. ■ 

Lemma 7. The expectation E[\T^ ] under the Levy fiight model 
is bounded by 

E[m] < 

LI U 
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Proof: Using Lemma |6] we can give a bound on E[[T]] in 
Lemma |4] as 

E[m] =f;P{T>r}<P„f;(Pr- 

r=0 r=0 

By (v) in Lemma |5] we have P < 1 for any r > 0. Thus, the 
expectation E[[r]] is bounded by the geometric series which 
converges to ■ 

Step 3 (Analysis of the delay scaling): In Lemma [1] we 
have analyzed the order of the per-node throughput A^r of 
the scheme tt. The results in Lemma ID (v) in Lemma |5] 
and Lemma |7] allow us to analyze the order of the average 
delay £>fr, which is shown in Lemma [8] 

Lemma 8. Let the communication range r scale as Q{n^) 
(0 < /3 < 1/4). Then, the average delay D^^ of the scheme tt 
under the Levy flight model with parameter a £ (0, 2] scales 
as follows: 

= 0(min(n(i+")/2-/3^^))^ 

Proof: Here, we provide a sketch of the proof with details 
given in Appendix B. Under the Levy flight model with 
parameter a € (0,2], we have (1 - P)-^ = e(n(i+")/2-/3) 
by (v) in Lemma |5] In addition, Pq = 0(1) by Lemma |3] 
Hence, from Lemma |4] and Lemma |2l we have 

< E[rT]] < = e(n(i+"'/2-'^). (23) 

In addition, under the Levy flight model, we have a trivial 
upper bound for all a £ (0, 2] as 

= Oin). (24) 

Combining (|23] | and (l24l i yields our lemma. ■ 

Step 4 (Analysis of the delay-capacity tradeoff): In the last 
step, we derive the delay-capacity tradeoff under the Levy 
flight model. By combining the capacity scaling in Lemma [T] 
and the delay scaling in Lemma [S] we get the following 
theorem. 

Theorem 1. Under the Levy flight model with parameter 
a G (0, 2], the delay-capacity tradeoff D{X) for per-node 
throughput A ~ &{n~^) (0 < rj < 1/2) is given by 

D{\) = 0(v/min(ni+"A,n2)). 

Proof: With the use of /3 = — log,, \/A, the scheme tt can 
achieve the per-node throughput Afr = A and the average delay 
Djj. = 0(-\/miii(ri^+"A, n'^)) by Lemma [T] and Lemma [8] 
respectively. Therefore, from (fTol ). we have our theorem. ■ 

VI. Delay Analysis for the i.i.d. Mobility Model 

In this section, we provide detailed analytical steps for 
obtaining the optimal delay under the i.i.d. mobility model. 
We again follow the four steps analogous to those used for 
the Levy flight model. 



Step 1 (Formulation of the average delay using the first 
meeting time distribution): From (fTTT i. the average delay 
under the scheme tt is obtained by 

Df, = P{d^I(s)}-E[min(r(,,d); I £ I{s)) \ d i I(s)] . (25) 

As pointed out in Step 1 of Section |V] the random variables 
T{i,d) for * G 2^(s) dependent. However, the dependency 
disappears when the nodes move to the next locations after a 
single time slot under the i.i.d. mobility model. The property 
of choosing a completely independent location at every time 
slot in the i.i.d. mobility enables this independence to occur 
By applying this observation, we derive a bound on r(i.c;) for 
i £ I(s) as follows: let |X(s)| denote the cardinality of the 
set X(s). We condition on the values of |I(.s)| and rewrite the 
expectation on the right-hand side of dZSl l as 

E[miii(T(,,,,i);*eI(s)) |d^I(s)] 

Tl-l 

= Y,V{\I{s)\^m\diI{s)} 

m— 1 

• E [ min (r(,.d) ; I e I(.s)) I |Z(s) I - m, d ^ I{s)] 
^P{|I(,s)| =m|d^I(s)}.E[min(T*,r*,...,T*)], 

m— 1 

(26) 

where T* {v = 1, . . . , m) denotes the first meeting time of 
the node d and the vth node in the set X(s), provided that 
|Z(s)| ~ m and d ^ I(s). Let Ti, . . . , T,„ be m independent 
copies of the generic random variable T. Then, we can derive 
a bound on T* in terms of as follows: 

T* ^ inf{t > I P(,;„,d)(t) < r, i„ e Z(s), d i X{s)\ 
< mi{t > 1 |L(,;„_d)(t) < r,i„ el{s),d(^l{s)} 

= l+T,. (27) 

Here, denotes the vth index in the set T(s) and = denotes 
"equal in distribution". The last equation comes from the 
aforementioned nature of the i.i.d. mobility model in which 
the locations of nodes are reshuffled at every time slot. 
We define a function f/ : {1, 2, . . . , n — 1} — > K by 

U{m) =E[min([T^,];u = l,...,m)]. (28) 

Note that discretization of a random variable to [T„] is 
for mathematical simplicity and it does not affect the result 
(i.e., order of the optimal delay) of this paper The function 
U{m) works as a tight upper bound on Z)# as shown in the 
following lemma. 

Lemma 9. The average delay D-a of the scheme n under the 
i.i.d. mobility model is bounded by 

Dir<Po + Po- E[U{B^„-2,Ps) + 1)] , (29) 

where Pq is defined in ( 1741 ). P^ = 1 — Pq, and P(„-2.P':) de- 
notes a binomial random variable with parameters n~2 (trial) 



and (probability). The function U (m) (m — 1, . . . , n — 1) 
can be obtained from the distribution of T by 

oo 

C7(m) = ^(P{r>r})™. 

r=0 

Proof Since T* < 1 + T^, < 1 + [T^,] for w = 1, . . . , m 
by (IZTT i. we have 

mill (T*, . . . , T„* ) < 1 + mill ( [Ti] , . . . , [T^l ) . 

By taking expectations, we have 

E[min(T*,...,r*)] < l + L/(m). (30) 

Since Xi{Q) is independent and identically distributed across 
node index i, each node i s) belongs to the set X(s) 
independently of each other with probability P^. Thus, the 
random variable |I(s)| — 1 (here, 1 is subtracted to exclude 
the case s € I(s)) subjected to the condition d ^ I(s) follows 
a binomial distribution with parameters 7i — 2 and P^, i.e., 

P{|2:(.s)| ^m\di I{s)} = P{B(„_2.P„c) = m - 1}. (31) 

By applying (l30l l and (l3Tl i to ( |26] |. we have 

E[min(T(,,d);zeI(s)) 

< 1+E[^(i?(„_2,pc) + 1)]. (32) 

Combining ^ and (O yields (|29] l. 

Since the random variable min([Tu] \v = 1, . . . , to) takes 
on only nonnegative integer values, U{m) can be obtained by 

oo 

t/M =^P{iiiiii(rT,l;w = l,...,m) >f}. (33) 

f =1 

By noting that [r„] is independent and identically distributed 

across u = 1, . . . , m, we have 

P{miii(rr„l;« = l,...,TO) >f} 

= (P{rrl>r})'" = (P{r>f-l})'", (34) 

where the second equality comes from the property that 
P{[T] > f } = P{r > f - 1} for all f = 1, 2, . . .. Hence, 
applying (|34] | to (|33] | and replacing f — 1 with r give the 
lemma. ■ 

Step 2 (Characterization of the first meeting time dis- 
tribution): In this step, similarly to the approach for the 
Levy flight model, we first analyze the characteristics of the 
function H{k, Iq) in (fT2] i under the i.i.d. mobility model (See 
Lemma [Toll. By exploiting the characteristics, we then derive a 
bound on the first meeting time distribution (See Lemma [TTTi. 
This bound enables us to derive a formula for the function 
[/(•) used in Lemma |9] (See Lemma [T2]|. 

As will be shown below, the characteristics of H(k, l^) 
under the i.i.d. mobility model are similar to those under the 
Levy flight model. Hence, an upper bound on P{T > r} can 
be derived using the probabilities P (= H{l,2y/Ti)) and Po 
also for the i.i.d. mobility model. The main difference is that 
the formula for H{1,Iq) is of different form and has a different 
scaUng property when = 2v^- 




Fig. 5. The illustration of S*{lo)- 

Lemma 10. Under the i.i.d. mobility model, the function 
H{k,lQ) in ([12\l has the following characteristics: 
( i) Let AX be a generic random variable for Xi (t) — Xj [t 



representing a location differential between nodes i and j 
Then, geometrically the function H(1,Iq) can be viewed as 
the probability of the location differential falling into a set 
S*{Iq) (c K^) defined for Iq £ {r,2y^] as 

S*{lo) = {a; eM^I ((0,/o),a;)nX>,(0) = 0}, (35) 

where the definitions of (•,•) and T>.(-) can be found in 
Lemma\5\ An example of S* (Iq) is shown in Fig. \5\ The set 
S*(lo) has a connection with the function H(l, Iq) as follows: 

H(l,l„)^P{AX eS*{lo)}. 
(ii) The function H(l,-) is nondecreasing. 
(Hi) From (ii), we have H(l,lo) < P for all Iq G {r,2^/n\. 

(iv) For k = 2,3, . . ., each function H(k, Iq) is also bounded 
above by P. Thus, for all k Cz N and Iq G (r, 2y/n\, we have 

H{k,lo) < P. 

(v) P is bounded above and below for all n £ N by 
1 



P < 1 - -sin^i 

TT 

P > 1 

2n TT 



2^^: 
2r 5 



\2./n) 



(36a) 
(36b) 



/n TT \2^Jn 
Proof: See Appendix C. ■ 

Similarly to Step 2 in Section |V] we derive a bound on 
P{r > t} in terms of P as follows: from (iv) in Lemma [TOl 
and Lemma |2] we have 

P{T > r} < (P)"Po, for r= 1,2,.... (37) 

Since P{T > t] = Po for r = 0, the bound in ^ also holds 
for r = 0. The above result is summarized in Lemma [TT] 

Lemma 11. Under the i.i.d. mobility model, the CCDF of the 
first meeting time T is bounded by 

P{T>t}<{P)^Po, for t = 0,1,..., 

where P and Pq are defined in M8^ and ([23, respectively. 

Proof: Combining Lemma [2] and (iv) in Lemma [TO] gives 
Lemma [TT] The detailed derivation was described earlier 
in ([37]). ■ 

"'The existence of the generic random variable for Xi (t) — Xj (i) is proven 
in Lemma [m in Appendix C. 
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Using Lemma (TT] we can give a bound on the function 
U{m) in Lemma |9] as 

OO CX3 

U{m) = ^(P{T > r})™ < ^((P)")^ (38) 

By (v) in Lemma [TOl we have (P)™ < 1 for any r > 0. Thus, 
U{m) is bounded by a convergent geometric series and we 
summarize the resuh in Lemma [T2] 

Lemma 12. The function U (m) (to = 1, . . . , n — 1) defined 
in (|2S| ii bounded under the i.i.d. mobility model by 

t7(m) < ^^"l" . 

1 - (P)™ 

Proof: Combining Lemma|9] (v) in Lemma [TOl and Lemma [TTI 
gives Lemma[T2l The detailed derivation was described earlier 
in ill. ■ 

The bound in Lemma [12] is essentially the same format 
with that of the slotted contact model under the i.i.d. mobility 
model. The only difference is that P additionally considers 
intermediate meetings. 

Step 3 (Analysis of the delay scaling): In this step, we analyze 
the order of the average delay under the i.i.d. mobility 
model. To efficiently handle the expectation E[C7(P(„_2,p=) + 
1)] in Lemma |9] we derive a bound on the expectation as 
follows: first, we rewrite E[t/(P(„_2,p^) + 1)] by conditioning 
on P(„_2,p-) as 

n-l 

E[C/(P(„-2,p„e) + l)] = ^ C/("j)-P{P(„-2,pc)=m-l}. (39) 

m— 1 

We then decompose ( |39] l into two terms as 

E[C/(P(„_2,pc) + l)] 

= ^ C7(m)-P{P(„_2,pc) =TO-1} 

m— 1 



= m - 1} 



<U{1) P{B(n-2,pc) !} + [/( [7^'!), (40) 

m— 1 

where 7 is a constant in (0, 1) and 7r^ implies the 7 fraction 
of the average number of nodes within the communication 
range of a source node. In ( l40l i. we used the property that 
U{m) is a nonincreasing function of m. Hence, by Lemma |9] 
and (|40] |. the average delay Dj^ of the scheme tt under the 
i.i.d. mobility model is bounded by: 

D^<Po + Po- U{1) ■ P{B(„_2,P.) < [7'^'] - 2} 

+ Po-t/(r7r2]). (41) 

The results in Lemmas [3] and [12] and (v) in Lemma [TO] 
allow us to analyze the order of the average delay Dfr, which 
is shown in Lemma [T3] 



Lemma 13. Let the communication range r scale as Q{n^) 
(0 < /3 < 1/4). Then, the average delay Dj^ of the scheme tt 
under the i.i.d. mobility model scales as follows: 

= O(^max(0,l/2-3/?))^ 

Proof: Here, we provide a sketch of the proof with details 
given in Appendix C. 
Order of Pq: By Lemma [3] 

Po = e(i). (42) 

Order ofU{l): By (v) in Lemma [Tol we have (1 - P)-i ^ 
Hence, combining (|42] | and Lemma [T2l vields 

Po 



uii)< 



1 -P 



(43) 



Order ofU{\jr^^): By ((42]i, we have (Po)^'''''^ = 6(1). In 
addition, by (v) in Lemma [TOl we have (1 — (P)r7'"'l)-i = 
0(„max(o,i/2-3/3))^ Hcncc, Lemma [H] givcs 

U{\^r^^) < ^^°^[^' ^ = e(„n— (0,1/2-3/3)) (44) 

^' ' - 1 - (P)M ^ ) \ > 

Order o/P{i?(„_2,pc) < [7r^] — 2}.- By using Chernoff's in- 
equality, for any fixed 7 <E (0, 1/3) and n > , we have 



P{S(n-2,P^) < 7^'} < CXp 

< 2 



37' 

1 / n — 2 — 8771 \ 2 



2 V 3n 
3n N 2 



0(n-2^) 



which results in 

P{P(„_2,p^) < \it'^ - 2} = 0(n-2^). (45) 
Combining (l4Tli-(l45Tl gives the lemma. ■ 

Step 4 (Analysis of the delay-capacity tradeoff): In the 

last step, we derive the delay-capacity tradeoff under the 
i.i.d. mobility model. By combining the capacity scaUng in 
Lemma [T] and the delay scaling in Lemma [13] we get the 
following theorem. 

Theorem 2. Under the i.i.d. mobility model, the delay- 
capacity tradeoff D{\) for per-node throughput X = 
e(n-'') (0 < ?y < 1/2) is given by 

D{X) = 0(v/max(l,nA3)). 

Proof: With the use of /3 — — log„ \/A, the scheme tt can 
achieve the per-node throughput A^r = A and the average 
delay Dji ~ 0(-\/max(l, nX-^)) by Lemma [T] and Lemma [T3] 
respectively. Therefore, from ([Tot . we have our theorem. ■ 

VII. Concluding Remarks 

In this paper, we developed a new analytical framework that 
substantially improves the realism in delay-capacity analysis 
by considering (i) Levy flight mobility, which is known to 
closely resemble human mobility patterns and (ii) contact 
opportunities in the middle of movements of nodes. Using 
our framework, we obtained the first delay-capacity tradeoff 
for Levy flight and derived a new tighter tradeoff for i.i.d. 
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Fig. 6. Examples of T) n X>r(a;): Area(X> n X>r(a;)) is nonincreasing as 
\x\ increases (Left). Wlien \x\ = n (e.g., x = (01, 0)), Area(T' n X>r(a3)) 
is minimized and a lower bound on Area(X>nX'r((\/"'i 0)) can be obtained 
in terms of 9 = 2 008"-*- ( 2^ )' given in i50\ (Right). 



mobility. For Levy flight, our analysis shows that the tradeoff 
holds 5(A) = 0(^min(ni+"A,n2)) for A = e(n-'') (0 < 
ij < 1/2) as shown in Figs.|2](a),|3](a), and [3](b). Our result is 
well aligned with the critical delay suggested in ifTsl . For i.i.d. 
mobility, our analysis provides D{\) = 0(-\/max(l, nX^)) as 
shown in Fig.|2](b). These tradeoffs are especially remarkable 
in both Levy flight and i.i.d. mobility for the constant per- 
node throughput (i.e., A = B(l)) as they demonstrate that 
the delay can be less than 9(n), which has been widely 
accepted for most mobility models. Our future work includes 
(i) an extension of our framework to analyze the delay-capacity 
tradeoff under Levy walk and (ii) another extension to capture 
correlated movement patterns among nodes. 

Appendix A 
Proof of Lemma[3] 

By the definition of Pq in (fl4] i. we have 

Po = P{i(.,,)(0) > r} = 1 - P{L(,,,)(0) < r}. (46) 

Let Fxi{o){') denote the CDF of Xi{0). Then, by conditioning 
on the values of Xi{0), the probability P{X(ij)(0) < r} in 
(l46T l can be rewritten as 

P{i(..,)(0)<r} 

= / P{i(.,,)(0) < r|X,(0) = x}dFx,io)ix) 

J-D 
J-D 

= f P{Xj(0)eX>.(a:)}dFx,(o)N, (47) 

where the last equality comes from the independence be- 
tween Xi{0) and Xj{0). Note that, since Xj{0) e T> with 
probability 1 and Xj{0) ^ Uniform(X>), the probability 
P{Xj{0) G Tfrix)} in the integral in ^} is given by 

where Area(5) denotes the area of a set 5 C M^. An example 
of T>C\'Dr{x) is shown in Fig.|6] From the figure, it is obvious 



that Area(I?nX'r(a;)) is nonincreasing as x approaches to the 
boundary of the disc "D. Hence, ( |48] l is bounded above by 

Area(X>n X>,-(0)) 



P{Xj(0) e Vr{x)) < 



^(V^)2 

Area(X>,.(0)) 



r 
n 



(49) 



In addition, it is bounded below by 

Area(X>nX>^((V^,0))) 



P{X,{0)e-Dr{x)}> 



^2 a 



> 



irn zTT n Ztt 
where 9 = 2cos~^(2^) (See Fig. |6ll. Since y/n > r, we 
have > 2cos^i(^) — Hence, the inequahty in dSOb is 
further bounded by 

P{X,iO)e-Dr{x)}>^. (51) 
By substituting (|49] l and ( BTT i into ( |47] |. we have 

^<P{L,,)(0)<.}<^, 
which, combined with (l46b . gives 

l--<Po<l-^. 

n Sn 
Appendix B 

Proofs of Lemmas for the Levy Flight Model 

Here, we give detailed proofs of Lemmas |5] and [8] which 
are used for analyzing the optimal delay under the Levy 
flight model. To prove Lemma |5] we need the following 
Lemmas [Il[l5l and [16] 

Lemma 14. For i ^ j and fc G N, let 

Av-(,.,)(fc)^v^,(fc)-y,(fc), 

where V^k) (representing the kth flight of a node ■) is defined 
in Then, under the Levy fiight model, AV(i has the 
following properties: 

(i) AV(i j)(fc) is independent of Xu{t) for all u = 1, . . . ,n 
and t G [0,k — 1]. 

(ii) AV(i j) (fc) is identically distributed across pair index {i, j) 
and slot index k. Hence, we use AV to denote a generic 
random variable for A Vj^ j) (fc). 

( Hi) For V € M^, let Z.v denote the angle at vertex enclosed 
by the line (0,t)) and the positive x-axis. Then, the angle 
ZA V is a uniform random variable on the interval (0, 27r] 
and is independent of the length |AV|. 

Proof: (i) For any u = 1, . . . ,n, X„(i) {Q < t < k ~ 1) 
under the Levy flight model is completely determined by 
Fu{k - 1) ^ {X„(0), y„(l), . . . , K(fc - 1)} (by the rela- 
tions ([T]i and (|2]i). Since Vi{k) is independent of Fu{k — 1), 
it is independent of By the same reason, Vj(fc) is 

independent of X„(t). Therefore, the difference Vi{k) — Vj [k) 
is independent of X^{t). 

(ii) Since each of the flight angle 9u{k) and the flight length 
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Fig. 7. An example of vectors satisfying the condition in (54): for a 
given (ni, Vj , v), the angle ZAV(i j) is determined by the angle ZVj. Since 
ZVj ~ Uniform[0, 27r], we have ZAV"{, j) ~ Uniform[0, 27r]. 

Zu{k) is independent and identically distributed across node 
index u and slot index k, the random variable Vu(fc) (= 
{Zu{k) cos 6 u (k), Zu{k) sin 9 u{k))) is also independent and 
identically distributed across u and k. Therefore, the difference 
Vi{k) — Vj{k) is identically distributed across pair index (i, j) 
and slot index k. However, it is not necessarily independent 
across while it is independent across k for a fixed 
(iii) To prove (iii), it suffices to show that for any u > 0, 

P{ZAy(,,,)(fc) < e I |AV(,^,)(fc)| ^v} = ^, (52) 

where < 6 < 2tt. In the following, we will prove ( |52] |. 

For simplicity, we omit the slot index k in V.{k) and 
A V(i (fc) in the rest of this proof. By conditioning on the 
values of (|V^;|, |), we can rewrite the probability on the 
left-hand side of ( |52] | as follows: 

P{^Ay(,.,.)<^?||AV(,,,.)H^;} 

= / p{zA%,) <0|(|y,|,|y,|,|A%,)|) = («.,i>„^>)} 

■P{m\,\V,\) = {v,,v,)\\AV^,,,)\^v}d{v,,v,). 

(53) 

For a fixed v > 0, consider an event {(IKIi 1^1) = {'^i^'^-'j)} 
such that 

P{(I'^".MV,I) = {v,,v,) I |AV-(,,,)| = v}>0. (54) 

An example satisfying ( |54] l is shown in Fig. [T] Under the con- 
dition {\V,\,\Vj\,\AV^,,j)\) = iv„Vj,v), the angle ZAV(,j) 
is determined by the angle ZVj as the figure shows. Since 
ZV-j - Uniform[0, 27r], we have ZAV(^,^j) - Uniform[0, 27r]. 
That is, for < 6* < 27r we have 

P{ZA%,) <e\i\V.\,\V,l\AV^,j)\) - {v,,v„v)} = 

Since the above equality holds for any {vi , Vj ) satisfying 
([54l i for a given v, the probability in (l53l) boils down to the 
following: 

p{zAy(,,)<e||Ay(,,,)| = «} 

= / P{(|V«|,|V,|) = (^'.,z^j)||AV(..,)| = z;}d(«„i;,) 

J{v,,Vj) 

9 

~ 2^' 



This completes the proof. ■ 

Lemma 15. Suppose fc G N and I € (r, 2^Jn\. Then, for any 
sets £(•) C [0, 2^/ri\ satisfying 

P{L{k -1)^1, L{t) G C{t), < < < fc - 1} > 0, (55) 
we have under the Levy flight model the following: 

P{I{k) = I L(fc - 1) = /, L{t) e C{t),0 <t<k-l} 
^P{AV eS{l)}. (56) 
The definitions of AV and S{1) can be found in Letmna\5\ 

Remark 1. Before proving the lemma, we give a remark. 
Lemma [75] implies that the future states of a meeting process 
under the Levy flight model depend only on the state at 
the beginning of the current slot, not on the sequence of 
events that preceded it. In addition, the conditional probability 
distribution of the future state described above is time homo- 
geneous (i.e., the probability in ( 1561 ) does not depend on the 
slot index k). This restricted time homogeneous memoryless 
property enables us to derive a bound on the first meeting 
time distribution as a geometric form ( See Lemma 

Proof: For notational simplicity, we let 

F{k-l) = {L{t)(^C{t),Q<t<k-l} (57) 
satisfying ( |55] |. For i ^ j and i > 0, let 

AX(,,,){t)^X.,{t)-X,{t). 

For simplicity, we omit in AX(i Then, by condi- 
tioning on the values of ZAX{k — 1), the left-hand side of 
(|56l l can be rewritten as 

P{/(fc) = 0\L{k - 1) = l,T{k - 1)} 

= / P{I{k)^0\ZAX{k-l)^9,L{k-l) = l,J^{k-l)} 
Jo 

(iFzAX{k~l)\L{k-l)=l,J^{k~l){(^)7 (58) 

where i^zAX(fc-i)|L(fe-i)=i,.F(fc-i) (•) denotes the CDF of the 
random variable ZAX{k — 1) conditioned that L{k — 1) = 1 
and T{k - 1). Since L{k - 1) = \AX{k - 1)|, the joint 
condition ZAX{k — 1) ^ 9 and L{k — 1) = Z is equivalent 
to AX{k - 1) = le^'^, where e^^ = (cos 6*, sin 6*). Hence, the 
probability in ( |58] l can be expressed as 

P{/(fc) = I ZAX{k -1)^9, L{k - 1) = /, F{k - 1)} 
= P{/(fc) = I AX{k - 1) = le^\F{k - 1)}. (59) 

The key idea of the proof is to use the following equality: for 

any k € N, I € (r, 2^/n], 9 G (0, 2tt], and T{k - 1), we have 

P{/(fc) - I AX{k - 1) = /e^^ T{k - 1)} 

(60) 

= P{AV eS{l)}. 
By substituting the combined result of ( l59l l and ( l60l ) into dSSl ). 
we have the lemma. 

In the following, we show ( l60l l. We first consider the event 
{/(fc) = 0}. By definition, the event {/(fc) = 0} occurs if 
and only if L{t) > r for all t G {k — l,k], equivalently. 
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AX{k-l) 




AX{k) 



Fig. 8. An illustration of a meeting event during the kth slot: in case of the 
blue line, L(fc — 1 + 5) > r for all 5 6 (0, 1], i.e., I{k) = 0. However, in case 
of the green line, there exist multiple 5 S (0, 1] such that L{k — 1 + 6) < r, 
i.e., I{k) = 1. 




Fig. 9. An example of S(l, 6): when 6» = |, S{1, 9) is identical to S{1). 

L{k - 1 + (5) > r for all (5 e (0, 1]. Since X,{k - 1 + 5) = 
(1 - S)X,{k - 1) + SX,{k) by ©, we have 

AX(fc -1 + 6) = {1- S)AX{k - 1) + SAX{k). (61) 

This implies that the event {L{k—1+S) > r for all 6 £ (0, 1]} 
occurs if and only if the following event occurs (See Fig. |8]i: 

|(AX(fc - l),AX{k))r\T>riO) = 0|. (62) 

We next consider the event {/(fc) = 0} conditioned by 
AX{k - 1) = le^^ and J"(fc - 1). Then, since X,{k) ^ 
X,{k - 1) + V^{k) by ©, we have 

AX{k) = AX{k - 1) + AV{k) = /e^^ + Ay(fc). 

Thus, given the conditions AX{k — 1) = le^^ and F{k — 1), 
( |62] i is reduced to the following: 

I (AX(fc - l),AX{k)) n X>,.(0) = 0| 
= |(ZeJS,?eJ» + Ay(fc)) nX»^(0) = 0} 
= {(0,AV(fc)) nX>^(-Ze^''') = 0| 
= {Ay(fc)e5(?,0)}, 

where 



S{1, d)^{x£B?\ (0,a;) n Vri-le^") = 0}. 

An example of S{l,d) is shown in Fig. |9] Hence, the proba- 
bility on the left-hand side of ( I6OI 1 becomes 

P{/(fc) = I AX{k~l) = ?e^^ J-(fc-l)} 
= P{AV{k) e S{l,e)\AX{k-l) = /e^^ J'(fc-l)}. (63) 



By (i) in Lemma [T4l Ay(fc) is independent of AX{k - 1) 
and F{k — 1), and thus we have 

V{AV{k) e S{1, 9) I AX{k - 1) = le'\F{k - 1)} 

= P{AV(fc) e 5(;,6')}. (64) 

In addition, by (ii) in Lemma [T4l 

?{AV{k) e 9)) = ?{AV e 9)]. (65) 

Finally, by (iii) in Lemma [141 the probability in ( |65] ) is invari- 
ant for any 9 G (0, 27r]. When 61 = f , we have S{1, f ) = S{1). 
Hence, the following holds for any 6* G (0, 27r]: 

P{Ay g5(/,6')} =P{AVg5(0}. (66) 

Combining ( |63] |. ( l64l i. ( |65] |. and (l66T l gives (|60] |. This com- 
pletes the proof. ■ 

Lemma 16. Lef Zi, Z2 and 9i , be independent copies of 
the generic random variables Z (flight length) and 9 (flight 
angle), respectively. Suppose that there exist constants c (> 0) 
and Zth (> 0) such that 



P{Z > z\ = — , for all z > z,i,. 



(67) 



Then, for all z > 2z,h we have 
ci 



where 



— <P{ZiCos9i - Z2COS6I2 > z} <—, 
z" 



ci^— (cos79)"flfi9(>0), 



27r 



2^+"c 



(cos7?)"c/i9(> 0). 



TT JO 

Proof: First, we will show that the distribution of Zcos9 is 
of the following power-law form; 

P{Zcos6' > X-} = — , for2:>Zth, (68) 



where ci = ^ J^^ {cos 1})" di) {> 0). By conditioning on 
the values of the random variable 9 ~ Uniform[0, 27r], the 
probability P{Z cos9 > x} can be rewritten as 

1 /-^^ 

P{Zcos9 > x} = — J P{Zcosi9 > x}di? 

= - / P{Zcosi9 > x}di?, (69) 

Jo 

where the second equality comes from the symmetry of the 
function cos with respect to = tt. For x > 0, the integral 
in ( |69l ) can be expressed as 



[ P{Zcosi9 > x}dz9 
Jo 



P{Zcosi9 > a;}di? + / P{Zcos^ > x}M, (70) 



where e G (0, -j). The first integral in dTOl i becomes 



P{Zcosd>x}M= ' P{Z > 

Jo cos iS 

= 4 (cos79)"d,?, (71) 

^ JO 
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where the first equaHty comes from cos i!) > for i? e [0, |- — e] 
and the second equaHty comes from ( |67l ) since -^^^-^ > Zth for 
X > Zfh. The second integral in ( |70] i is bounded by 

0< I P{Zcosd > x}M < I ldi? = e. (72) 

Combining (|69ll, (|70l), (|2B, and gives 

c f^^'' 

/ (cos t?)" di? < P{Z cos 6* > x} 



< 



Letting e ^ on (l73l l yields 



(cos^?)"dz9+ -. (73) 



P{Zcos6l > x} ^ / (cosz?)"dz?. 
Hence, we have 

P{Zcos0 > X-} — — , for X > Zth, 



where Ci = /^^ (cos??)" d?9 (> 0). This proves 

In the following, we derive the distribution of the random 
variable Zi cos Oi — Z2 cos by using ( l68l l. Since the event 
{ZiC0s6'i < §} n {Z2 cos 6*2 > — 1} implies the event 
{Zi cos — Z2 cos 6*2 < ^}, we have 
P{Zi cos 01 - Z2 cos 6*2 > z} 



z z 

< P{Zicos6'i > 2 ^2 cos 6*2 < 

< 2P{Zcos6' > 



(74) 



where the last inequality comes from the union bound and 
the symmetry of ZcosO (i.e., ZcosO = —ZcosO). Suppose 
z > 2zth. Then, by applying (l68T l to ( f74] i. we further have 

2"+ici c„ 



P{Zi cos 6*1 - Z2 cos 6*2 > z} < 



(75) 



Similarly, since the event {Zi cos 6*1 > z} n {Z2 cos6'2 < 0} 
implies the event {Zi 0036*1 — Z2 0036*2 > z}, we have 
P{Zi cos 9i - Z2 cos 02 > z} 

> P{Zi cos 01 > z and Z2 cos 02 < 0} 
= P{Zi cos 01 > z}P{Z2 cos 02 < 0} 

= ^ = -, (76) 

2z« z" 

where the first equality comes from the independence between 
Zi cos 01 and Z2 cos 02, and the second equality comes from 
(l68T l and the symmetry of Z2 cos 02. Combining (iTST i and (|76] | 
gives the lemma. ■ 

Proof of Lemma|5] 

A. Proof of (i) 

By choosing k — I, I — Iq, and £(0) = (r, 2-^71] in 
Lemma [TSl we have 

P{/(1) = I L(0) = /o, i(0) e (r, 2^^} = P{AV e 5(Zo)}. 
Since {L(0) = Iq} n {i(0) € {r,2y/^]} = {L(0) = ^n}, we 
further have 

P{/(1) = I L(0) = lo} = P{AV e S{lo)}. (77) 
By definition, H{1, Iq) = P{/(1) = | L{0) = Iq}. Thus, we 
have from (|77]i that H{1, lo) = P{AV £ S{lo)}. 



B. Proof of (ii) 

Suppose r < lo < h < ^^/n. Then, it is obvious from the 
definition of S{-) in ^ that S{lo) C 5(^1) (See Fig. |4|. 
Hence, we have 

P{AFe5(;o)} <P{AFe5(;i)}, 

which is equivalent to lo) < ^/^(l, ^1) by (i) in Lemma|5] 

C. Proof of ( Hi) 

By (ii) in Lemma|5] we have H{1, lo) < H{1, 2^) (= P) 
for any G (?', 2^\. 

D. Proof of (iv) 

Recall the definition of H{k, lo) for = 2, 3, . . .: 

H{k, lo) ^ p{m = o|/(fc- 1) = . . . = /(I) = 0, m = ^o}. 

By conditioning on the values of L{k — 1), the probabihty 
H{k, lo) can be rewritten as 

HikJo) 

= / P{/(fc)=0|L(fc-l) = ?,/(fc-l) = ... = /(l)=0, 

Jr+ 

L{0) = ?o}dFi(fc_i)|/(fe„i)=.. -/(i)=o,L(o)=io(Oi (78) 

where FL{k_i^\j(^k-i}=...=i{i)=o,Lio)=io{-) denotes the CDF 
of L{k - 1) conditioned that /(fc - 1) = . . . = /(I) = and 
L{0) = /g- Here, we integrate I = L{k — 1) over (?■, 2^/71] 
due to the condition I{k — 1) = 0. By using Lemma [TSl the 
probability in the integral in ( iTST i is simplified as follows: 

P{/(fc) = 0|L(fc- l) = Z,/(fc- l) = ... = /(l) = 0,L(0) = /o} 
= P{/(fc)=0|L(fc-l) = ?,L(t)e (^,2^^], 

< < < fc- l,i(0) = ^o} 
= P{AV e S{1)}. (79) 

By (i) and (iii) in Lemma |5] the probabihty P{AV G S{1)} 
is bounded for all I E (r, 2^^] by 

P{Ay G 5(/)} = i7(l,0 < A (80) 

By substituting the combined result of ( |79] l and dSOb into ( fTSl ). 
we have for all fc = 2, 3, . . . and lo E {r, 2^/n\ the following: 

HikJo) <P-l = P. 
This proves (iv) in Lemma |5] 

E. Proof of (v) 

By (i) in Lemma H P = P{AV G S{2y^)}. To derive a 
lower and upper bound on P, we define a subset «S^(2y^) 
and a superset •S^(2y^) of the set S{2^/n) as depicted in 
Fig. [To] Then, we have 

P{AV e 5^(201)} < P < P{AF e 5+(2V^)}. (81) 
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Fig. 10. The geometric definitions of the superset S'^{2y/n) (Left) and 
the subset S~ {2yfn) (Right) of the set <S(2y'?i) used in the proof of (v) in 
LeiTima|5] 



By (iii) in Lemma [TH the probabilities P{Ay G 5=^(2^)} 
are obtained by (double sings in same order) 

P{AV e S^{2^)} 

= 1 -P{Ay i 5±(2V^)} 

= l-P{|Ay|>2V^±r}.^i|^, (82) 

where Lp{2y/n) is the central angle associated with S^{2^/n) 
(See Fig. [Toll. From the geometry in Fig.[TO] the angle Lp{2y/n) 
is given by 



(f{2\/n) = 2 sin 



2v^ 



(83) 



We now consider the probabilities P{|AV| > 2^/n ± r} in 
(l82T i. For notational simplicity, we denote AV = (A14, AVy). 
Then, by (ii) in Lemma [T4l 

AV^ = Z^{k) cos6',(fc) - Zj{k) cos 9 J (k), 
AVy = Z,{k) sin 6,{k) - Zj{k) sin6'j(fc). 

Note that for any v = {vx,Vy) e and r/ > 0, \vx\ > i] 
implies \v\ > i], and \v\ > rj implies > r]/V2 or \vy\ > 
r]/\/2. Hence, P{\v\ > rj} is bounded by 

P{|^^..| >ri}< P{\v\ > V} (84a) 
< HK\ > v/V2 or \vy\ > v/V2}- (84b) 

Since 9i{k) and Oj{k) are independent and uniformly dis- 
tributed over (0, 27r], AVx is symmetric, i.e., AVx = —AVx- 
Thus, applying ( I84al i with v = AV and r/ = 2^^ + r yields 

P{|Ay| > 2^^: + r} > P{|A14| > 2y/^ + r} 
= 2P{AVx > 2VH + r}. 

Since Zth in Lemma [16] is a constant independent of n, there 
exists a constant nth,; G N such that 2^^ + r > 2zth for all 
n > ?^th,;- Hence, by Lemma [T6l we have for all n > nth,; 



P{|Ay| > 2^/^ + r} > 2ci 



1 



2J^+r 



(85) 



Since cos 61 = sin6i for 9 - Uniform[0, 27r], |AT4| ^ \AVy\. 
Thus, applying ( |84bt with ?; = A V and jy = 2y^ — r yields 

P{|AV| > 2V^-r} 

< P{|AF^| > {2y^ - r)/V2 or \AVy\ > {2^/li - r)/V2} 

< 2P{\AVx\ > (2^/H-r)/V2} 
= 4P{A14 > (2V^-r)/^/2}. 

By the same reason as above, there exists a constant nth,ti G N 
such that (2y^ — r)/\/2 > 2zth for all n > nth.u- Hence, by 
Lemma [T6] we have for all n > nth.u 

/ 1 



P{|Ay| >2V^-r}<2"/2+2c ^ ^ 

Combining ([SB, ([82li, ([83ll, ([Mil, and ^ yields 

2q / 1 \" t / r 

P > 1 



(86) 



TT \2-\/ri + r 

2"/2+2c„ / 1 



2v^y ' 



TT \2y^ — r 
for all n > nth = niax(nth,;, nth,u)- 



2y/?i 



Proof of Lemma[8] 

To complete the proof of Lemma [Sj it remains to show that 
(i) (l-Py^ = e(n(i+")/2-^) and (ii) Df, = 0{n). Without 
loss of generality, we assume r = (0 < /3 < 1/4). 

A. Proof of (i) 

To prove (i), we need the following: for any x £ [0, 1], 



X < sm (x) < —X. 



(87) 



The proof of Wt\ is given at the end of this section. From 
( |19a| i in Lemma [5] with r ~ n^, we have for all n > rith the 
following; 



2q 

1 - P > — 



1 



TT V2\/n + n^ 



/3-1/2 



Since ii-y— e [0,1] for any (3 e [0,1/4] and n e N, we 
further have from the lower inequality in ( [87] ) that 

ci ( 1 



1 -P > 



TT \2y^ + n 



,/3-l/2 



Hence, we have 



(l-P)-i 2"7r 
limsup , ,„ „ < < cxj, 

which gives 

(1 - P)-i = (9(n(i+")/2-'5). (88) 

Using a similar approach as above, from ( |19b| ) in Lemma [5] 
and the upper inequality in ( [87b . we have for all n > nth the 
following: 

1 



1 - P < 2"/2c„ 
Hence, we have 



2Jn-nl^ 



„(1+q)/2-/3 

limsup ^ < 2 "'^Cu < oo, 
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which gives 

(l-P)-i = (89) 
Combining dSSl l and ( [89] l proves (i). 

Proof of W7\: For \x\ < 1, the function sin~^ (x) can be 
calculated using the following infinite series; 

oo 
l=Q 

where di = 41 i^p^2''^2i+i) Hence, for any x £ [0, 1], we 

have a lower bound on sin^^ (x) as 

sin"^ (x) >dox^x. (90) 

Since x'^'^"^ < x for all / = 0, 1, . . . and x G [0, 1], we have 

00 

sin^^ (x) < X di. 

Note that J^Zo = sin"^ (1) = f ■ Hence, for any x G [0, 1], 
we have an upper bound on siii^^ (x) as 

sin-i (x) < ^x. (91) 

Combining ( |90l l and ( |9T] ) proves (Isfl l. 

B. Proof of (ii) 

Without loss of generality, we assume P{Zq, > Zth} = 1- 
(In this proof, subscript a is added to all random variables to 
specify the underlying parameter a of the Levy flight model.) 
Then, from dSll, we have P{Zq > z} ^ (^)" for all z > Zth, 
which gives for any < ai < a2 < 2 and z > z± the 
following; 

nZc, >z}= (^)" > (^)" =P{Z^., > z}. (92) 
The inequality in ( |92] i shows that for any t2 > ti > having 
a sufficiently small difference e ^ t2 — ti > 0, we get 

P{iai(i2) > r-|i„,(ti) > r} < P{L^,{t2) > r\L^,{h) > r}, 
which results in 

P{Tc., >t2\T^,>ti}< P{T„, > t2 I > ti}. (93) 
Note that since P{T„, > t} ^ P{T„, > t,Ta^ > t - e} for 
t > e, we can express P{Ta^ > t} in a nested form as 

P{Ta, >t}= P{T^, >t\T^,>t- (}P{T^, >t-e}. 
Using the nested form continuously, we have 

P{r„, >t}= P{T^, >t\T^,>t-e} 

X P{T„, > t - e I >t~2e} 
X . . . 

X P{T,, >t - \ t/e\e\T^, > 0} 
X P{T„, > 0}. (94) 
Hence, by applying ( |93] l to ( |94l i, we have 

PlTai > t} < P{T„, > 1 1 > t - e} 

X P{T„, > t - e I Ta, > t - 2e} 
X . . . 

xP{T„, >t - [t/e\e\T^, > 0} 

X P{T^, > 0}. (95) 



Note that P{Tqi > 0} = P{Lai(0) > r}. In addition, since 
X,;(0) ^ Uniform(X>) for alH = 1, . . . , n regardless of a, we 
have P{T„, > 0} = PjTa^ > 0}. Thus, the right-hand side 
of (|95] l boils down to PlT^^ > t}, and consequently 

P{T„, >t}< P{r„, > i} for all t > 0. (96) 

Due to the property in (|96] l, the average delay under the Levy 
flight model with a parameter a E (0, 2) is dominated by the 
one under Brownian motion {a = 2), which is shown to be 
0{n) Ei, i.e., 

D^ = 0{n) for an a e (0,2]. 



Appendix C 

Proofs of Lemmas for the i.i.d. Mobility Model 

Here, we give detailed proofs of Lemmas [TO] and [T3] which 
are used for analyzing the optimal delay under the i.i.d. 
mobility model. To prove Lemma [TOl we need the following 
Lemmas [HI [IS and [19] 

Lemma 17. For i =^ j and t > 0, let 

AX(,,,)(t)^X,(t)-X,W, 

where X. (t) denotes the location of a node ■ at time t. Then, 
under the i.i.d. mobility model, AX(j [t) has the following 
properties: 

(i) AJC(i (k G N) is independent of Xu{t) for all u = 

and t (z [0, fc — 1]. 

( ii) AJC(-j j-j (t) is identically distributed across pair index {i, j) 
and time t (> 0). Hence, we use /S.X to denote a generic 
random variable for AJf(i 

(Hi) The angle ZAJC is a uniform random variable on the 
interval (0,27r] and is independent of the length \/S.X\. 

Proof (i) For any u = l,...,n, (0 < t < fc - 1) 

under the i.i.d. mobility model is completely determined by 
Gu{k - 1) ^ {X„(0), . . . , Xu{k - 1)} (by the relation O). 
Since Xi{k) is independent of G„(A; — 1), it is independent of 
Xu{t). By the same reason, Xj{k) is independent of Xu{t). 
Therefore, the difference Xi{k) — Xj{k) is independent of 

Xu{t). 

(ii) For any i 7^ j and t > {), Xi{t) and Xj{t) are inde- 
pendent and identically distributed. Therefore, the difference 
Xi{t) — Xj{t) is identically distributed across pair index 

and time t. However, it is not necessarily independent neither 
across nor across t. 

(iii) To prove (iii), it suffices to show that for any x > 0, 

P{ZAX <e\\AX\=x} (97) 

where < 6* < 27r. By noting that ZXi{t) Uniform[0, 27r] 
for any i — 1, . . . ,ti and t > and using a similar approach 
as in the proof of (iii) in Lemma [14] we can prove (iii) in 
Lemma [TT] Due to similarities, we omit the details. ■ 
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Fig. 11. 



An example of S*{1,6): when 9 = ^, S*{1,0) is identical to 



Lemma 18. Suppose fc e N and I G (r, 2^yn\. Then, for any 
sets £(•) C [0, satisfying 

P{L{k -l) = l, L{t) e C{t), < t < fc - 1} > 0, 

we have under the i.i.d. mobility model the following: 

P{I{k) = I L(fc - 1) = /, L{t) e C{t), 0<t<k^l} 

= P{AX e S*{1)}. (98) 

r/ie definitions of AX and S*{1) can be found in Lemma [7q1 

Remark 2. Before proving the lemma, we give a remark. 
As Lemma \T5\for the Levy flight model. Lemma [TSl implies 
that the future states of a meeting process under the i.i.d. 
mobility model depend only on the state at the beginning of 
the current slot, not on the sequence of events that preceded 
it. In addition, the conditional probability distribution of the 
future state described above is time homogeneous (i.e., the 
probability in (|9S| does not depend on the slot index k). This 
restricted time homogeneous memoryless property enables us 
to derive a bound on the flrst meeting time distribution as a 
geometric form ( See Lemma [77]). 

Proof: Using a similar approach as in the proof of Lemma [TSl 
we can prove Lemma [18] The difference is that the key idea 
of this proof is to use the following equality: for any fc e N, 

I e (r, e (0, 27r], and T{k - 1), we have 

P{Iik) = 0\AX{k - 1) = le^^,T{k^l)} 

= P{AXg5*(0}, (99) 

where the definition of J'{k — 1) can be found in ( fSTb . Then, 
similarly to the proof of Lemma [15] using the key equality in 
i99[ we can prove Lemma [18] Due to similarities, we omit 
the details. 

In the following, we show ( i99] l. We first consider the event 
{/(fc) = 0}. Since (i6Tl i also holds for the i.i.d. mobility model, 
by the same reason in the proof of Lemma [15] the event 
{/(fc) = 0} occurs if and only if the following event occurs: 



AX{k - 1), AX(fc)) nX»^(0) = 0|. 



(100) 



We next consider the event {/(fc) = 0} conditioned by 
AX{k — 1) = le^^ and T{k — 1). Under these conditions. 



dlOOI l is reduced to the following: 

I (AX(fc - 1), AX(fc)) n -Drio) 

= |(ZeJ», AX(fc)) nX»^(0 
= {AX(fc)e5*a,^)}, 



) = 0} 



where 



S*il,e) ^{xeM.'^l {lei^,x)r]T>riO) = 0}. 

An example of S*{l,d) is shown in Fig. [TT] Hence, the 
probability on the left-hand side of (i99l ) becomes 

P{J(fc) = 0| AX(fc-l) = le^^^Tik^l)} 
= P{AX{k)eS*{l,e)\AX{k-l) = le^^,T{k-l)}. 

(101) 

By (i) in Lemma [Tt! AX{k) is independent of AX{k — 1) 
and F{k — 1), and thus we have 

P{AX(fc) e d) I AX{k - 1) = le^\ F{k - 1)} 

= P{AX(fc) e 5*(/,e')}. (102) 

In addition, by (ii) in Lemma [TT] 

P{AX(fc) e S*{1, 9)} = P{AX e S*{1, e)}. (103) 

Finally, by (iii) in Lemma [17] the probability in (11031 ) is 
invariant for any G (0,27r]. When 6* = -I, we have 
S*{1,^) = S*{1). Hence, the following holds for any e 

(0,27r]: 

P{AX eS* (1,0)} = P{ AX eS*{l)}. (104) 

Combining ([ToTT i. ( [T02] |. ([T03] l. and ( [T04] | gives ([99]). This 
completes the proof. ■ 

Lemma 19. Let /|ax|(') denote the probability density func- 
tion of the random variable \AX\ under the i.i.d. mobility 
model. Then, it is bounded by 
2x 

f\AX\{x) < — , for all X e[Q,2y/n]. 
n 

Proof: We will prove this lemma by showing the following: 
P{x - § < |AX| < a;+ f 



lim sup ■ 

£^0 



(105) 



From (ii) in Lemma [17] we have |AX| = |AX(,_j)(0)|. 
Hence, by conditioning on the values of Xi{0), the probability 
in ( 11051 ) can be rewritten as 



Pix - -< \AX\ < X 



■<|AX 



1} 

,-)(0)|<x+||x,(0) = M}dFx,(o)(t 



/ P{X, (0) e (u) I X, (0) = It} dFx, (0) (u) 

J-D 



where 72. 



{x,e} 



< \v — u\ < X - 



s}- By 



independence between Xi{0) and Xj{0), we further have 



{.-1<|AX|<.+ 1} 
= / P{X,(0)e7^(,,,)(«)}d^^x.(o)N- 

J-D 



(106) 
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Note that, since Xj (0) G X> with probabiUty 1 and Xj{0) ^ 
Uniform(X>), the probability P{Xj{0) e 7?.(^. ,)(«)} in the 
integral in ( |106b is given by 

Area(X>n72.(^,,)(M)) 



P{x,(0) = 



< 



Area(X>) 

Area(7?.(^^,)(M)) 



(107) 



In addition, for any u E T> and sufficiently small e (> 0), the 
area Area(7?.(2. e)(tt)) is calculated as 

Area(7^(,,,)H) = <j \ ' 

if X ~ 0, 

(108) 



^7r(x + f)2, 
27ra;e, if a; > 0, 



By applying the combined result of (1107b and (IIO8I 1 to (1106b 
we have 



p{:r-|<|AX|<x+|}<<j2 
which gives 



if a; > 0, 



if a; = 0, 



P{a;-f <|AX|<a:+f} ^ If, if a; > 0, 
lim sup — 



< 

eio t 1^0, if a; = 0, 

_ 2x 
n 

This proves the lemma. ■ 
Proof of Lemma[To] 

A. Proof of (i) 

Similarly to the proof of (i) in Lemma |5] we can prove (i) 
in Lemma [To] by applying Lemma [Ts] with fc = 1, / = Iq, and 
£(0) = {r,2y/n\. Due to similarities, we omit the details. 

B. Proof of (ii) 

Suppose r < Iq < li < Then, it is obvious from the 

definition of 5*(-) in ^ that 5*(Zo) C S*{h) (See Fig. E). 
Hence, we have 

P{AX e 5*ao)} < P{AX e <S*(/i)}, 

which is equivalent to Zq) < ^^(1: '1) by (i) in Lemma [TOl 

C. Proof of ( Hi) 

By (ii) in Lemma[Tol we have H{1, Iq) < H{1, 2^) P) 
for any Iq S (r, 2\/n]. 

D. Proof of (iv) 

By following the approach in the proof of (iv) in Lemma |5] 
we can prove (iv) in Lemma [TOl based on Lemma [TS] and (i) 
and (iii) in Lemma [TO] Due to similarities, we omit the details. 



y 






^ — tA ^ 


1 





Fig. 12. The geometric definition of it is tlie central angle of the arc 

{x 6 5*(2Y^)||a;| = x\ depicted in red. 



E. Proof of (v) 

By (i) in Lemma [TO] and (iii) in Lemma [iTl we have 
P = P{AX e S''{2^/^)} 



27r 



-/|AX|(a;)da;, 



(109) 



where 0(a:) is the central angle of the arc {x^S* {2y/n)\\x\ = 
a;} (See Fig. fTZb. and ./|ax|(') is defined in Lemma [T9l From 
the geometry in Fig. [12] the angle 4){x) is given by 

^(x) = 2cos-(^)+2cos- Q 

= 27r-2sin"^ f^^) -2sin"^ f-V (110) 
\2y/n/ \x/ 

where the second equality comes from the identity cos^^{d) = 
f -sin-^(6') (-f < 61 < f ). By substituting ( ITTOI l into ( [T09b . 
we have 

1 



P = I 1 sin" 

TT 



2^ 



sm 

TT J,.+ \X 



2^^ 

1 



P{|AX| > r} 



f\AX\{x)dx. 



(Ill) 



Based on (II 1 lb . we derive an upper bound on P as follows: 



P < [1 sin , ^ 

- V TT ^2^/n 

1 1 / r 
< 1 sin"^ 



P{|AX| > r} 



TT \2^.y 

This proves the upper bound in ( I36ab 



Using (111 lb again, we derive a lower bound on P as follows: 
since |AX| = |AX(0)| by (ii) in Lemma [TT] and |AX(0)| = 
L(0) by definition, we have P{|AX| > r} = P{L(0) > r}. 
Hence, by Lemma [3] the probabiUty P{|AX| > r} in (lUlb 
is bounded by 



P{|AX|>r} = Po>l-^. 



(112) 



By Lemma [19] the integral in (111 lb is bounded by 

/ sin-M-j/|Ax|(a;)da: < - / 



sin I — |a;dx. 

a'/ 



(113) 
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Let y = r/x. By the change of variables, the integral on the Order of U{\^r'^~\): To complete the proof of (l44l i. it remains 

right-hand side of iU3\ is solved as to show (1 - {py^'"^^ = 0(„max(o,i/2-3^s)-) p^j. jj^j^^ 

2 r'^Vn , ^ , will show the followings; 

n X+ IxJ'^'^'' (iii) (1 - (P)hr-^)-^ = o(„max(0,l/2-3«)^ 

= + ^./i_l!+4sin-i f^V (114) (i^) (1 - (^')^^'-'^)-' = 17(n--(°-V2-3«)_ 

2n -v/n V 477. \2-i/ri/ , 



From ( 1115b . we have 



By applying (11121) . ( I113l l. and (11141) to (llUb . we have 

^2\pn) I \ n ) To simplify ( II 17b . we will use the following bound: for any 



P> l--sin"M — ^ -ll ^ 27r 



TT 



^2 2r / 4 _i / r \ a; G [0, 1] and y > 0, 
sin 



2n ttV^V 4n tt V2Vn/ 1 - x^^^^ = (1 - + x + . . . + xr^^T"^) 

> 1 - ^ - ^ - ^ sin- (-^). > (1 - 



This proves the lower bound in (I36bb . 

Proof of LemmaHJ] 
Order ofU{l): To complete the proof of (l43T l. it remains to 



>{l-x)yxy. (118) 

By applying (1118b with a; = 1 - "'^l'^ (e [0, 1]) and y ^ 
7r^ = 771-^'^ (> 0) to the right-hand side of dl 17b . we have 

i-(p)r7'-^i > fi-!^ y . 



show (1 - P)-i = 6(771/2-^). For this, we will show the Hence, we have 

followings: (l - (P)^-^'-'!)"^ 2tt / n/3-i/2s -^n^" 



limsup TPTTf, < — limsup 1 

(i) (1 - = 0(77^/^"^), n^/^-iP 7 „_^oo V 27r / 

(ii) (1 - Py = n{n^/'^-^) 



(119) 

— 1 / 2 2(3 

To obtain the order of (1 — ^-^ ) , we take a logarithm 

Without loss of generality, we assume r = 77/^ (0 < /3 < 1/4) ^^^^^-^^ ^^^^^ ^^^^^^^ ^j^^ ^-^^^^ behavior: 

in the rest of this appendix. From ( I36ab in Lemma [Tol with ^i-i/2 -jn^f 

r = n^, we have for all 77 g N the following: lim log I 1 



2tt 

l-^>--n^(^). =limlogl-" 



27r 

Since e [0,1] for any /3 G [0,1/4] and 77 € N, we fur- / „,3-i/2 , „i/2-/3 



7 lim ( ri^/^^^/^iog ^ 

n^oo \ \ 



ther have from the lower inequahty in dSTi t (i.e., x < sin ^ (x) n^oo \^ V 27r 

for any X e [0, 1]) that 3 2 , nP-^/^Y^'"' 

= -7 hm 77-''^ 'l^- hm logll 

71' n~>oo \ ZTT / 

1-^^^- (115) . f lim 773^-V^. (120) 

, 2n n— i-oo 

Hence, we have /3-1/2 23 

- Hence, we have lim„^oo log(l — ^^^^ — )-7" — g for f3 E 



limsup — ^ ■„ < 27r < 00, [0, 1/6). That is. 



n^oo 



^/3-l/2^„^„2,3 



which proves (i) (1 - P)-i = 0(771/2-''). J*™, V"^ 9^J " ^^^^^ 

Using a similar approach as above, from (ElHi in Lemma[lO] By combining (HUll and dnB, for /3 e [0, 1 /6) we obtain 

and the upper inequahty in dSTl i (i.e., sin^^ (x) < -ja; for any / ^p\\l■r'^^\^^ 2 

a; e [0, 1]), we have for all 77 G N the following: limsup -^^ \ ,i ' — < — < 00, 

- n^(^-^ 277/^-1/2 5 , /77'3-i/2x ,., / ' 

I- p < - ^ — \- _ sin"i ( - ) which results in 

2J-V2 "5^7/3-1/2 ^ (l-(P)rT'-'T)"' =0(771/2-3^), for /3e [0,1/6). (122) 



< 



77 



2 ' TT ' 4 From (Iinil, we have 1 - (P)rT''-'l > 1 - (1 - ii^2F-)'i'"'^ 

^ -|- 7^,(3-1/2 (116) Hence, we have 

h ' (i-(p)r^^^i)-' 

Hence, we have lim V y_J_ [ — 

1/2 — 3 o -7 n— >oo 77 

' '^2 7 1 ('17^^ 

hmsup < - + 7 < 00, ^ 1 ^^^•'^ 



(l-P)-l " TT 4 



< 



7n 

lim 1 _ f 1 - ' 

which proves (ii) (1 - = ^{7^^-^) umn^oo ( 



27r 



2fi ■ 
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The Umit lini„^oc(l - "^^^^V'" is obtained from ( fT20] i as 
follows: 

lim log 1 = - lim log 1 — 

n— >-oo \ ZTT / n— )-oo \ zTT 



That is. 



lim ( 



^ 27r ) 



\~oo, if /3 e (1/6, 1/4]. 

fexp(-2^), if/3 = l/6, 
10, 



if /3e (1/6,1/4]. 

(124) 

By substituting ( 11241 ) into ( 11231 ). for /? G [1/6, 1/4] we obtain 
lim sup ■ _ 

if /3e (1/6,1/4], 

< oo, 

which results in 

(1 - (P)^'''^'!)"^ = (9(n°), for /3 € [1/6,1/4]. (125) 



Combining ( I122l i and ( |125t proves (iii). 
From dl 16l l, we have 



(126) 



.TT 4/ 

To simplify ( 1126b . we will use the following bound: for any 

X G [0, 1] and y > 0, 

i-x^'^T = {\ - x)(\ ^ X ^ . . . ^ x^y^-^) 

<{\-x)\y\. (127) 
By applying WTT\ with a; = 1 - (| + Dn/^-^/^ [g, 1]) and 
y ~ = 771^'' (> 0) to the right-hand side of ( |126l l. we 
have 

Vtt 4/ 

Hence, we have 

„l/2-3/3 2 7\ 

lim sup : =- < ( h - [7I < 00, 

n->oo (i_ (P)r7r^l)-i - Vtt 4/'" 

which results in 

(1 - (P)r7-^1)"' = f7(ni/2-3/3)^ for ^ g [0^ 1/4]. (128) 

In addition, since we have 

limsup r < 1 < oo, 

n^oo (l _ {p)\lr^^y^ ~ 

which results in 

(1- (P)^'"-'^)"' = r!(n"), for /3e [0,1/4]. (129) 
Combining ( I128I I and ( |129t proves (iv). 

Ort/er of P{B^^n~2.p^) < [t^^I ^ 2}.- By Chernoff's inequal- 
ity, the lower tail of the distribution function of the binomial 
random variable i?(„_2.p<:) for x < [n — 2)P^ is bounded by 

P{i3(„_,,p,) < 4 < exp ( - ^^'^[''^"pf ) ■ (130) 

From Lemma [3J we have = 1 — Po > Suppose < 
7 < I and n > (or, equivalently, n — 2 > S-fn). Then, 



we have (n — 2)Po > (n — 2)|^ > 3"fn^ = jr"^. Hence, we 
can apply x = 7r^ to ( |130l l under the conditions < 7 < ^ 



and n > ., ^ , and we obtain 

— 1— 37' 



((n- 2)P„=-7r2)2 



P{P(„_,P,)<7.^}<exp^- 

Since > ^ and n — 2 > 8771, the term {n — 2)P^ — jr 
in ( 11311 ) is bounded below by 



(131) 

2 



(n-2)P„^-7r^ > (n-2) — -7r 



3n 



2 n — 2 — 3771 
3n 



(> 0), 



from which we have 

{in-2)P^-^r^)- > 



2\ 2 fn — 2 — 37n\ 2 ^ 



(132) 



From Lemma [3] we also have P^ = 1 — Po < 7^- Hence, the 
term 2(77, - 2)P„= in ([mi l is bounded above by 2(n - 2)P„= < 
2(71 — 2)— < 2r^, from which we have 



1 1 

2(71 - 2)PS ~ 2^' 



(133) 



Thus, by (I132l i and (1133b . the argument of the exponential 
function in ( 1131b is bounded below by 

((77-2)P„^-7r2)' ^ 1/ 71-2-37^ x2 ^ 
2(71 - 2)P„^ - 2 V 371 / ^ ' 
which gives an upper bound on ( 1131b as 

1 / 71 — 2 — 3711 \ 2 



P{B(„-2.p„=) < ir'} < exp (^-- ^ 
Since exp(— a:) < i for all x > 0, we further have 



P{Bin-2,Ps) < ir^} < 2 
Therefore, for r = 71*^ we have 



371 



71 — 2 — 3771 



lim sup ?T7^ < 
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< CXD, 



71-2/3 - (1 _ 3^)2 

which results in P{P(„_2.pc) < 77-^} = 0{n~'^^). Since 
P{B(„„2,p„^) < \ir^} - 2} < P{P(„_2.p^) < ir^}, we have 

P{S(„-2,p„^) < \lA - 2} = 0(71-2^*). 
This completes the proof. 
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